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ABSTRACT 


Ongoing  work  seeks  to  apply  the  technology  of  resonance  ionization  mass  spectrometry 
(RIMS)  to  problems  related  to  nuclear  forensics  and,  in  particular,  to  the  analysis  and 
quantification  of  the  debris  from  nuclear  detonations.  As  part  of  this  effort,  modeling  and 
simulation  methods  are  being  applied  to  analyze  and  predict  the  potential  for  ionization 
by  laser  excitation  of  isotopes  of  both  uranium  and  plutonium. 

Early  work  focused  on  the  ionization  potential  of  isotopes  of  uranium,  and  the 
present  effort  has  expanded  and  extended  the  previous  work  by  identifying  and 
integrating  new  data  for  plutonium  isotopes.  In  addition  to  extending  the  effort  to  this 
important  new  element,  the  work  described  in  this  thesis  implemented  more  accurate 
descriptions  of  the  spatial  distribution  of  the  laser  beams  to  improve  the  accuracy  of 
model  predictions  compared  with  experimental  results  as  well  as  an  ability  to  readily 
incorporate  new  experimental  data  as  they  become  available. 

The  model  is  used  to  estimate  ionization  cross  sections  and  to  compare  the 
relative  excitation  potential  for  two  isotopes  as  a  function  of  wavelength,  irradiance,  and 
bandwidth. 
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I.  INTRODUCTION 


A  discovery  that  happened  over  Christmas  vacation  in  Vienna,  1938,  has  had 
global  impacts  on  the  world  of  nuclear  physics  [1].  Physicists  Lise  Meitner  and  Otto 
Fritsch  were  contemplating  a  finding  by  Otto  Hahn  and  came  across  a  phenomenon  that 
was  previously  thought  impossible:  the  splitting  of  a  uranium  nucleus.  Fission,  as  Frisch 
called  it,  was  bom.1  Soon  after  this  discovery — credited  also  to  Fritz  Strassman — 
scientists  realized  these  reactions  also  emitted  secondary  neutrons.  These  secondary 
neutrons  could  in  turn  enable  chain  reactions,  releasing  vast  amounts  of  energy  (see 
Figure  1).  This  energy  could  be  harnessed  and  used  for  a  virtually  endless  supply  of 
clean,  electrical  power.  With  the  world  at  war,  however,  more  of  the  focus  was  directed 
toward  the  alternate  use  of  the  energy  release:  the  destructive  capabilities  of  nuclear 
fission  and  atomic  weapons. 


Figure  1.  Uranium-235  Fission  chain  reaction,  from  [1]. 


1  Enrico  Fermi  had  detected  fission  in  1934,  when  he  began  bombarding  elements  with  neutrons,  vice 
the  more  popular  protons,  and  found  the  production  of  new  elements.  However,  he  was  not  fully  aware  of 
the  implications  until  after  the  discovery  of  the  German  scientists  [1], 
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A.  THE  ATOMIC  BOMB 

The  availability  of  such  vast  amounts  of  energy  during  a  period  when  the  world 
was  rapidly  approaching  global  conflict  resulted  in  international  scientists’  urgent  efforts 
to  harness  this  power  and  use  it  to  arm  their  respective  nations  with  a  deadly  weapon. 
During  the  first  few  years  of  World  War  II,  Nazi  Germany  was  thought  to  be  in  the  lead 
in  the  nuclear  arms  race.  They  had  discovered  fission  in  1934  and  from  the  spring  of 
1940,  had  dedicated  a  large  part  of  the  Kaiser  Wilhelm  Institute  in  Berlin  to  uranium 
research  [1]. 

In  the  United  States,  in  August  of  1939,  a  letter  written  by  physicist  Leo  Szilard  in 
consultation  with  Edward  Teller  and  Eugene  Wigner,  and  signed  by  Albert  Einstein,  was 
sent  to  President  Franklin  D.  Roosevelt  [1].  The  letter  discussed  the  threat  that  Gennany 
might  develop  atomic  weapons  and  urged  that  the  United  States  begin  its  own  nuclear 
weapons  effort  [1].  The  letter  ultimately  resulted  in  the  fonnation  of  the  National  Defense 
Research  Committee  (NDRC)  and  the  Manhattan  Project,  in  which  the  scientific 
community  conducted  the  research  in  the  United  States  in  their  own  quest  for  an  atomic 
weapon. 

At  the  time,  Ernest  Lawrence,  a  Professor  at  the  University  of  California  and  a 
key  player  in  the  Manhattan  Project,  had  been  involved  in  experiments  using  a  cyclotron, 
as  shown  in  Figure  2.  Research  involving  neptunium  revealed  that  the  element  decayed 
into  another  trans-uranium  element  [1].  This  element  was  to  be  known  as  plutonium.  This 
discovery  was  crucial  because  it  suggested  that  there  was  a  possibility  of  producing  large 
amounts  of  fissionable  plutonium  from  neutron  irradiation  of  uranium-238,  and 
chemically  separating  it.  The  discovery  of  this  phenomenon  and  the  use  of  plutonium  in 
modem  nuclear  weapons  are  important  factors  of  direct  relevance  to  this  thesis  research. 


2 


Figure  2.  Ernest  Lawrence  prepares  his  cyclotron  in  the  radiation  lab  at  UC- 

Berkeley,  from  [2]. 

In  April  1945,  President  Roosevelt  passed  away,  and  Vice  President  Harry 
Truman  ascended  to  the  presidency.  Truman  was  not  involved  in  any  of  the  secrecies 
associated  with  the  Manhattan  Project  and  therefore  had  to  be  thoroughly  briefed.  In  the 
war,  Germany  was  on  the  brink  of  surrender  and  Japan  was  being  decimated.  However,  it 
was  a  general  consensus  that  Japan  would  fight  to  the  very  end.  Some  of  the  American 
policy  makers  thought  the  only  way  to  put  an  end  to  the  fighting  was  to  deliver  the 
crushing  blow  of  an  atomic  attack  on  recognized  strategic  cities. 
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Figure  3.  Trinity  Bomb,  0.053  seconds  after  detonation,  from  [1]. 


Around  the  same  time,  the  Americans  were  preparing  to  test  a  plutonium  weapon 
in  addition  to  the  uranium  bomb  that  was  well  developed.  The  Trinity  bomb  was 
detonated  on  July  16,  1945  [1].  This  test  meant  the  United  States  could  now  use  a  second 
type  of  bomb  against  the  Japanese.  This  was  true  start  of  the  Atomic  Age. 

B.  NUCLEAR  FORENSICS 

Nuclear  forensics  is  neither  a  secret  technology  nor  a  new  one.  Records  show  it 
dates  back  to  1949,  when  samples  drawn  by  the  United  States  from  high  altitudes  off  the 
coast  of  China  proved  that  Russia  had  detonated  a  nuclear  device  in  Central  China  [3]. 
However,  in  more  recent  times,  the  nuclear  smuggling  and  trafficking  phenomenon  has 
led  to  the  development  and  improvement  of  nuclear  forensics  [3]. 

Nuclear  forensics  is  a  fairly  new  science,  whose  modem  birth  event  is  considered 
to  be  the  seizure  of  Find- 1  in  Augsburg,  Gennany,  consisting  of  72  contraband  uranium 
pellets  (see  Figure  4)  [4].  Standard  techniques  used  to  analyze  these  pellets  included: 
potentiometric  titration  for  detennining  the  uranium  content,  thermal  ionization  mass 
spectrometry  (TIMS)  for  measuring  isotropic  composition  and  optical  microscopy  for 
determining  the  macroscopic  parameters  [4]. 
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Figure  4.  Uranium  pellet  from  Find  1,  from  [4], 

If  the  United  States  were  attacked  by  a  nuclear  weapon,  questions  would  arise  [5]: 

•  What  is  the  material  in  the  weapon? 

•  How  bad  is  the  damage?  And  how  much  worse  will  it  get? 

•  Who  did  it?  And  did  they  have  help? 

•  Are  there  more  out  there? 

•  What  should  I  do  about  it? 

These  questions  must  be  addressed  in  the  context  of  an  enormous  amount  of  stress 
and  pressure  on  the  President  to  give  the  public  accurate  answers  in  a  short  amount  of 
time.  These  answers  will  provide  the  President  with  information  required  to  make 
appropriate  and  decisive  decisions. 

In  the  case  of  plutonium,  we  know  it  can  be  released  from  various  sources,  such 
as  nuclear  weapon  detonations,  nuclear  processing  facilities  or  reactor  accidents  [6]. 
Therefore,  isotope  ratios  can  reveal  important  information  in  regards  to  a  weapon’s 
origins.  Unfortunately,  the  plutonium  concentration  can  be  low  in  environmental  samples 
of  debris,  so  sensitive  testing  methods  are  necessary.  Resonance  ionization  mass 
spectrometry  (RIMS)  has  been  proposed  as  an  alternative  or  supplement  to  existing 
analytical  methods  to  improve  our  nuclear  forensics  capabilities  [6]  and  will  be  described 
further  in  Chapter  II.  An  important  feature  of  this  technology  is  its  ability  to  detect  ultra¬ 
trace  amounts  of  plutonium.  As  seen  in  Table  1,  the  ratios  detennined  by  RIMS  compare 
well  with  previously  published  values. 
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Table  1.  Isotope  ratios  of  plutonium  in  uranium  ammunition,  from  [6], 


Fraction  of  Pit  as  measured  by  Fraction  in  weapon- grade  plutonium 
RIMS  [%]  (3  a  error)  [%]  according  to  Fetter  et  al.  [4] 


Pu-239 

93.6(5.9) 

93.5 

Pu-240 

6. 2(0. 7) 

6.0 

Pu-241 

0.09(0.03) 

0.44 

Pu-242 

0.06(0.01) 

0.015 
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C.  OVERVIEW  OF  THESIS  RESEARCH 

Ongoing  work  seeks  to  apply  RIMS  technology  to  problems  related  to  nuclear 
forensics  and,  in  particular,  to  the  analysis  and  quantification  of  debris  from  nuclear 
detonations.  As  part  of  this  effort,  modeling  and  simulation  methods  are  being  applied  to 
analyze  and  predict  the  potential  for  ionization  by  laser  excitation  of  both  uranium  and 
plutonium  isotopes.  Early  work  focused  on  the  ionization  potential  of  isotopes  of 
uranium,  and  the  present  effort  has  expanded  and  extended  the  previous  work  by 
identifying  and  integrating  new  data  for  plutonium  isotopes.  The  present  work  has 
implemented  a  more  accurate  description  of  the  spatial  distribution  of  the  laser  beams. 
This  will  improve  the  accuracy  of  model  predictions  when  compared  with  experiment 
results.  It  will  also  assist  in  the  ability  to  readily  incorporate  new  experimental  data  as 
they  become  available.  The  model  is  used  to  estimate  ionization  cross  sections  and  to 
compare  relative  excitation  on  two  isotopes  as  a  function  of  wavelength. 

Detennining  laser  power  for  optimal  ionization  is  crucial  for  detennining  the 
spatial  distribution  of  each  of  the  excited  states  of  the  elements  of  interest.  The  impact  of 
small  shifts  in  ionization  peaks  can  be  overcome  by  the  use  of  broader  band  lasers.  This 
will  also  enable  simultaneous  and  unifonn  ionization  of  multiple  elements.  Careful  and 
detailed  understanding  of  irradiance  levels  and  predicted  isotopic  responses  is  required. 
This  research  effort  has  sought  to  define  these  levels  and  explore  ways  to  further  refine 
these  results. 

Previous  efforts  were  conducted  to  model  the  laser  ionization  potential  of  isotopes 
of  both  uranium  and  plutonium  for  potential  nuclear  forensic  applications  that  involve  the 
use  of  RIMS  to  quantify  debris  from  nuclear  detonations.  The  research  presented  here 
represents  an  expansion  on  that  previous  work  by  identifying  and  integrating  new  data  for 
plutonium  isotopes  into  our  modeling  package  with  the  intent  of  gaining  a  more  accurate 
representation. 
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D.  OUTLINE 

Chapter  II  introduces  RIMS  and  the  physics  involved  in  the  process.  It  also  briefly 
describes  isotope  shifting  and  some  of  the  causes.  Chapter  III  introduces  the  rate  equation 
model.  It  describes  in  detail  the  components  of  the  model.  This  is  use  to  further 
understand  the  relative  ionization  probability  of  Pu.  It  also  allows  for  model  testing 
without  using  resources  which  are  used  experimentally.  Chapter  IV  presents  our  model 
results  and  compares/contrast  them  with  experimental  data.  Chapter  V  is  the  summary  of 
the  modeling  comparisons.  It  also  provides  direction  for  future  research. 
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II.  RESONANCE  IONIZATION  MASS  SPECTROMETRY 


Resonance  ionization  mass  spectrometry  (RIMS)  uses  the  energy  of  photons  from 
lasers  to  excite  and  ultimately  remove  an  electron  from  the  ground  state  of  an  atom  [7]. 
This  process  starts  by  creating  a  cloud  of  material  from  a  debris  specimen  by  laser 
desorption.  This  cloud  of  material  is  then  subjected  to  excitation  by  multiple  pulsed  lasers 
finely  tuned  to  excite  a  transition  from  one  state  to  the  next  for  the  specific  elements  of 
interest.  The  target  atoms  are  then  promoted  above  their  ionization  threshold  and  are 
extracted  and  sent  through  a  time-of-flight  mass  spectrometer  and  counted  [7]. 

As  seen  in  a  simplified  example  in  Figure  5,  element  A  has  two  states:  an  excited 
state  with  non-zero  and  finite  lifetime  (given  by  t)  energy,  and  a  ground  state.  The 
probability  for  finding  the  excited  state  at  FWHM  is  given  by  T,\.  If  the  wavelength  (/.)  is 
close  to  this  energy  difference,  then  the  probability  of  excitation  and  de-excitation  is 
high.  Isselhardt  explains,  “If  the  wavelength  varies  greatly  from  the  energy  difference,  as 
seen  for  element  B,  then  there  is  a  low  probability  of  excitation’’  [8]. 


Figure  5.  Energy  level  diagrams  of  two  random  two-state  atoms  of  different 

elements,  from  [8]. 
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A. 


THE  PHYSICS 


The  following  provides  an  outline  of  the  physics  involved  in  RIMS.  Isselhardt,  in 
[8],  provided  the  foundation  for  the  topic.  The  following  is  a  summation  of  his  work  in 
both  of  these  references. 

1.  Atomic  States 

Two  atomic  states  (Figure  6)  can  be  characterized  by  the  difference  in  their 
energies,  E2  —  Et  —  hcj0  >0.  If  an  atomic  state  can  decay,  it  must  therefore  be  time- 
dependent.  The  time  dependent  Schrodinger  equation  can  be  used  to  derive  the  time 
dependent  wave  equation  T  (t).  Using  the  normalized  wave  function  v|/0,  which  is 
stationary  in  the  time  independent  Schrodinger  equation,  we  can  show  from  [8] 

_.E 

'f'(t)  =  iyQ  .  e  lht 

In  its  upper  state,  the  atom  may  decay  with  an  average  lifetime  x,  where  x  =  7  . 

A 


Figure  6.  Two-state  atomic  model,  from  [8]. 


If  an  atomic  state  can  decay,  then  the  probability  of  it  occurring  in  T  is  given  in  [8]  by 

(W\W)  =  e~u 
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This  requires  the  wave  function  to  be  complex.  The  energy  for  an  atomic  state  that  can 
decay  is  defined  as  E  —  E0  +  iy,  -2 y  and  must  equal  h/z,  which  will  be  referred  to  as  V. 
This  represents  the  FWHM  of  the  energy  distribution  about  Eg.  When  it  is  expressed  as  a 
function  of  angular  momentum,  T  becomes  1/x,  since  E  —  hen.  The  time  dependent 
Schrodinger  equation  can  now  be  written  as 


'f'(t)  =  W0  ■  e  h 


The  Fourier  transform  allows  us  to  describe  a  function  of  time  in  terms  of 
frequency.  Also,  “  E  =  hen  and  dE  —  hden,  the  probability  for  finding  a  time-dependent 
state  in  energy  is 


1 


r 


with  P  (E)  representing  a  Lorentzian  distribution  with  a  FWHM  of  f,”  shown  in  Figure  7. 
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Figure  7.  Lorentzian  Function  with  f  =  0.2,  from  [8]. 
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2.  Emission  and  Absorption 

The  line  shapes  of  spontaneous  decay  and  induced  emission  and  absorption  are 
almost  identical.  Their  probability  distributions  are  closely  related.  The  electromagnetic 
radiation  density  in  a  cavity  can  be  described  by  Planck’s  distribution  and  given  in  [8]  as, 


u(w )  dw 


u(v) 

2n 


dv  — 


hw3  dw\ 

n2c3  J 


1 


Isselhardt  describes  these  variables  as  “it(w)  dw  is  the  energy  per  unit  volume,  h 
is  the  reduced  fonn  of  Planck’s  constant,  c  is  the  speed  of  light,  k  is  Boltzmann’s 
constant,  and  T  is  temperature  in  Kelvin”  [8].  This  equation  is  legitimate  when  a  balance 
is  reached  between  energy  absorbed  at  a  given  frequency  by  atomic  vapor  and  the  energy 
emitted  by  the  vapor  at  the  same  frequency  range.  These  two  populations  are  related  by 


N2 

Ah 


—  hWg 

e~kT~ 


The  probability  that  an  atom  will  either  emit  a  photon  into  or  absorb  a  photon  from  a 
radiation  field  by  decay  or  emission  can  be  given  as 

rl  P 
Llr  abs 


dt 
d  P 

ur  emit 

dt 


=  B12  •  uw(w) 
=  B2 1  •  uw(w) 


dP, 


decay 


dt 


=  A 


21 


According  to  Isselhardt  in,  “uw(w)  is  the  energy  density  of  the  radiation  field  as  a 
function  of  angular  frequency.  Ajj,  Bi2,  and  B2i  are  Einstein  coefficients  for  decay 

emission  and  absorption”  [8].  Since  state  2  can  only  decay  into  state  one,  A21  must  equal 

1 

the  inverse  of  the  lifetime  of  state  2  (—  =  r2).  Using  this  knowledge  combined  with  the 

t2 

fact  that  at  thennal  equilibrium,  these  quantities  can  be  expressed  as  [8] 

(A2 1  +  B21  ■  uw(w))  •  N2  =  B12  ■  uw(w)  •  Ah 
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Substituting  Planck’s  distribution  into  the  equation  yields  [8] 


aw(w)  = 


A21  1 


Comparing  this  result  to  Planck’s  distribution,  the  relationship  for  Einstein’s  coefficients 
can  be  derived  as  follows  [8]: 


02^12  —  01^2 1 
hw 


A 


21 


n2c3 


B 


21 


a.  Angular  Momentum  and  Degeneracy  States 

Energy  degeneracies  can  be  removed  in  the  presence  of  an  electric  or  magnetic 
field.  However,  in  thennal  equilibrium,  the  total  population  of  a  given  level  will  relate  to 
the  number  of  degenerate  states  with  a  level.  Using  a  two  state  model,  this  can  make 
things  extremely  difficult,  if  we  were  to  consider  the  degeneracy  states  of  both  these 
levels.  All  atoms  have  vast  quantities  of  levels  to  be  considered;  therefore  accounting  for 
the  transitions  between  these  levels  must  be  considered.  If  we  were  to  consider  these 
transitions  as  incoherent,  we  can  apply  a  rate  equation  approximation.  In  heavy  saturation 
conditions,  from  [8],  we  can  use 

N 2  _  02 
Ni  9i 

Atoms  at  each  state  will  now  be  considered  equal  to  their  degeneracy  levels,  which  will 
be  important  when  a  resonance  ionization  scheme  has  to  be  considered. 


b.  Selection  Rules 

Conservation  of  momentum  and  parity  imply  there  must  be  rules  for  each  electron 
transition.  The  transition  dipole  matrix  is  stated  in  [8]  and  represented  by 

\Dik\e  =  e  •  ef  W^rWidq 

When  this  is  equal  to  zero,  it  is  said  the  transition  is  dipole  forbidden.  Since  the  parity  of 
the  dipole  operator  is  odd,  the  parity  of  the  wave  functions  involved  in  the  transition  must 
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be  odd.  Therefore,  the  rule  of  parity  came  to  be  known  as  Aparity  =  yes.  Other 
selection  rules  for  dipole  transitions  are  Am.  and  A] .  Their  values  are  listed  in  Table  2. 


Table  2.  Selection  Rules  for  Dipole  Transitions,  from  [8]. 


Component 

Rule 

Am, 

o,±i 

A] 

0,  ±1(0  to  0  forbidden) 

Aparity 

yes 

3.  Lineshape:  Doppler  Broadening 

Since  the  atoms  are  moving  in  various  directions  in  space,  we  have  to  consider 
they  have  a  frequency  vg  and  have  a  velocity  v.  These  atoms  move  in  a  parallel  direction 
of  propagation  with  the  laser.  A  Doppler  shift  can  be  expected  and  is  given  by  [8]  as 

vD  =  v0(l  +  ^j 

Thermalized  atom  velocities  in  the  gas  phase  can  be  described  by  the  Maxwell- 
Boltzmann  probability  distribution  given  by  [8]  as 

I  hr  -mv2 

PMdv=JMTen^dv 


Frequency  distribution  about  v„  is  given  by  [8]  as 


1  -Q-^o)2 

$(v)=-=e  2<A 

V  Ztccj 

and  is  known  as  a  Gaussian  distribution  with  a  standard  deviation  given  by  [8]  as 


This  gives  the  FWHM,  or  a  Doppler  width  of  given  in  [8]  as 


Avd  =  2  •  v0 


\2kT  ■  In  2 
mC 2 
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Using  the  constant  values  and  converting  to  wavelength,  the  following  expression  is 
derived  in  8  as 

AAd  =  7.16  x  1CT7  •  A0  • 

T  is  in  [K]  and  M  is  in  amu. 

4.  Ionization 

Resonance  excitation  from  the  ground  state  to  the  final  ionization  usually  requires 
several  intennediate  pathways.  Possible  transitions  are  pictured  in  Figure  8.  Lifetimes  for 
the  first  excitation  steps  last  approximately  10ns.  The  cross-section  for  absorbing  a 
photon  is  given  by  o  —  A2 /2n,  which  is  only  for  an  atom  at  rest.  The  excitation 
probability  is  determined  using  dW(t )  =  er/(t)  dt,  with /(t)  *  101S  photons/  (cm2s)  and 
corresponds  directly  with  a  power  of  100mW/cm2  [6].  These  values  easily  represent 
saturation  of  the  first  transition.  However,  according  to  Wendt  in  [6],  these  values 
represent  an  ideal  situation  of  constant  velocities  in  all  atoms;  an  effective  cross-section 
must  be  used.  This  is  the  reason  direct  non-resonant  ionization  (pathway  of  the  left  side 
of  Figure  8)  is  less  than  favorable.  The  cross  section  provided  would  be  approximately 
10’17cm2.  This  process  is  less  than  favorable  for  actinides  with  complex  atomic  spectra 
(plutonium)  and  can  be  avoided  by  heavily  populating  the  auto-ionization  state,  which  is 
a  low-lying  state  above  the  ionization  potential.  This  state  converts  energy  into  an 
electron  and  photon  and  decays  within  Ins.  For  excitation  into  the  Rydberg  state, 
ionization  probability  is  greatly  increased  [9]. 
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Rydberg 

States 


Ground  State 


Figure  8.  Generic  Ionization  Scheme  for  a  multi-step  RIMS,  from  [6], 


Figure  9  is  a  partial  energy  diagram  of  plutonium  showing  the  ionization  scheme 
of  Pu  and  Pu.  It  includes  the  energy  of  the  levels  and  the  wavelengths  of  the  lasers 
used  to  excite  the  atoms  to  the  next  state.  This  scheme  shows  plutonium  atoms  excited 
from  the  ground  state  (E  =  Ocm'1)  to  their  first  excited  state,  with  energy  levels  of 
23765.75  cm'1  and  23766.4  cm'1  [9].  These  atoms  are  then  excited  to  the  second  state,  to 
energies  of  35568.85  cm'1  and  35568.85  cm’1  [9],  Finally,  they  reach  ionization  through  a 
relatively  narrow  (6.03eV)  auto-ionization  state  with  an  energy  of  48597.65  cm'1  [9]. 


3-color  Pu-RIMS  (with  isotope  .shifts) 


Figure  9.  Plutonium  3-color,  3-photon  ionization  scheme,  from  [9]. 
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5. 


Autoionization 


Autoionization  was  first  identified  by  Ugo  Fano  in  1961  where  he  describes  states 
involving  two  excited  atoms  with  Ejot  >  ionization  potential  [10].  Atoms  in  the 
autoionization  state  can  either  decay  radiatively  to  less  than  the  ionization  potential  or 
through  non-radiative  decay  by  forming  an  electron  and  ion.  Isselhardt  represents  this 
form  of  decay  by  i(n1n2)  — >  f(n3k )  and  describes  the  variables  as  “nj  is  the  bound  state 
of  the  electron  and  k  is  the  state  of  the  free  electron.  i(n1n2)  is  a  discrete  bound  state  with 
Ejot  >  ionization  limit  ”  [8].  Partially,  in  [10]  the  width  depends  on  the  coupling  of  states 

rai  =  2n\{i\V\f)\2 

and  is  known  as  the  partial  width  of  decay  to  the  continuum.  V  is  the  interaction  matrix 
between  the  first  and  final  states.  Figure  10  gives  a  generic  autoionization  process. 

Doubly 


Figure  10.  Autoionization  schematic,  from  [8]. 
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Autoionization  resonant  lineshape  (Fano  profile)  [10]  is  described  by 


1  +  £z 

2  (E-E0) 

£  =  - 

r 

Dl 

D2R12 

Fano  defines  the  variables  as  “era  is  the  absorption  cross  section  for  forming  a  bound 
state,  ob  is  the  cross  section  for  ionization  into  the  continuum,  £  is  the  reduced  energy  for 
the  doubly  excited  state,  D{  is  the  amplitude  to  the  doubly  excited  state,  Dz  is  the 
amplitude  to  the  continuum  state,  and  R12  is  the  transition  rate  between  these  states”  [10]. 

B.  ISOTOPE  FACTORS 

Shifting  of  isotopes  is  caused  by  many  factors.  Two  of  these  will  be  explored. 

The  nuclear  mass  difference  factor  usually  involves  the  lighter  isotopes. 
According  to  Isselhardt,  it  results  from  the  fact  that  “the  additional  mass  of  a  neutron  is  a 
large  part  of  the  total  mass.  This  will  lead  to  a  small  difference  in  the  mass  of  a  bound 
electron,  which  will  in  turn  affect  the  energy  of  the  excited  state”  [8]. 

The  volume  difference  factor  is  due  to  the  electrical  charge  distribution  for 
different  isotopes  and  affects  the  heavier  isotopes.  Since  the  atomic  mass  difference  is 
relatively  small  compared  to  the  total  atomic  mass,  with  each  additional  neutron,  the 
nuclear  radius  may  see  large  changes  in  its  nuclear  shape.  This  will  subsequently  lead  to  a 
change  in  the  electric  field. 

1.  Odd  Isotopes 

Isselhardt  explains,  “Odd  massed  nuclei  have  one  unpaired  nucleon.  Therefore, 
these  nuclei  have  an  angular  momentum  (/)  >0. 1  nuclear  couples  with  the  /  of  the  electrons 
in  its  orbit”  [8].  The  angular  momenta  of  the  nuclear  (7)  and  the  orbital  electron  (J)  must 
be  coupled  into  total  angular  momenta  (F). 
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The  left  diagram  (a)  in  Figure  1 1  shows  an  even  isotope  with  I  =  0,  s  substrate 
level  rrij  =  0  and  substrate  in  the  upper  level  mj  =  -1,  0,  1.  The  right  diagram  (b),  shows 
an  odd  isotope  with  the  nuclear  spin  angular  momenta  of  /  =  3/2.  In  contrast  with  the 
even  isotope,  there  is  two  states,  F  =  3/2  and  1/2  where  the  energy  is  shifted  by  the 
isotope  shift  and  split  by  the  Casmir  fonnula  given  in  [1 1]  and  seen  in  Figure  11. 


A 

AE(F )  =  —C  +  B 


3  C(C  +  1)  —  47/(7  +  1)(/  +  1) 
87/(2/  —  1)(2/  -  1) 


C  =  F(F  +  1)  —  /(/  +  1)  -J(J  +  1) 


a)  Even  Isotope 

1  =  0 


b) 


Odd  Isotope 
1=  1/2 
F  =  I  +  J 


J  =  I 


Airij  =  0 


nij  -1  0  1 


F  =  3/2 
F  =  1/2 


F  =  1/2 


-3/2  -1/2  1/2  3/2 


Figure  11.  Partial  level  diagram  for  even  and  odd  isotopes,  from  [8], 


Substituting  rrif  for  mj  and  taking  note  Am,j  =  0,  we  can  easily  see  there  are  four 
allowable  transitions.  Angular  momenta  differences  affect  the  population  rate  of  each 
substrate.  This  will  generate  differences  in  the  ionization  probabilities  for  both  isotopes. 
Using  AJ  —  + 1  avoids  fractionation  effects  of  changes  in  the  selection  rules  [12]. 

2.  Even-Odd  Effect 

Peter  Lambropoulos  and  Andreas  Lyras  were  modeling  the  even-odd  isotope 
ration  for  tin  (Sn)  as  a  function  of  excitation  laser  power,  when  they  noticed  differences 
in  the  ionization  probabilities  of  the  even  and  odd  isotopes  [13].  However,  if  the 
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transitions  are  saturated,  then  no  differences  should  be  noticed.  They  accounted  for  these 
differences  as: 


•  More  states  in  the  odd  isotope — the  number  of  states  that  can  be  excited 
by  changing  the  selection  rules  (F  instead  of  J).  This  will  lead  to 
fractionation,  increasing  the  number  of  atoms  that  are  present  during  the 
time  of  transition  saturation. 

•  Dilution  of  the  oscillator  strength. 

According  to  Payne,  Allman,  and  Parks  [12],  we  can  overcome  the  effects  of 
fractionation  by: 

1 .  Using  fast  excitation  lasers 

2.  Changing  the  bandwidth  of  the  lasers 

Both  of  these  methods  increase  the  likelihood  of  non-resonant  ionization  of 
unwanted  isotopes.  However,  if  the  transitions  are  sufficiently  saturated,  these  effects  can 
be  minimized.  Since  care  must  be  taken,  the  cross  section  of  the  ionizing  laser  is  kept 
larger  than  the  excitation  lasers,  providing  the  assurance  that  ions  will  be  fonned  outside 
the  areas  where  the  ionization  process  is  saturated. 

C.  FRACTIONATION  IN  IONIZING  PROCESS 

The  main  causes  for  the  discrepancies  in  ionization  probability  of  isotope  ratios 
are  described  below: 

1 .  Bandwidth — Isselhardt  explains,  “Ionization  probabilities  for  isotopes  will 
differ  depending  on  laser  irradiance  over  the  resonant  cross  section”  [8]. 

2.  Selection  Rules — Angular  momentum  differences  can  lead  to  a  divergence 
in  the  number  of  states  the  isotope  can  transition  to  [14]. 

3.  Dynamic  Effects — Ionization  not  saturated  can  lead  to  fractionation 
because  the  odd  isotopes  in  the  excited  states  when  saturate  are  greater 
than  their  even  counterparts.  Lambropoulos  and  Lyras  in  [13]  suggest, 
“Ionization  rates  may  differ  due  to  different  angular  momentum  quantum 
numbers.” 
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III.  RATE  EQUATION  MODEL 


The  model  that  is  applied  in  this  research  was  originally  developed  to  authenticate 
experimental  data  resulting  from  tests  performed  on  isotope  ratios  to  modifications  in 
laser  bandwidth.  It  was  used  to  calculate  ionization  probabilities  of  uranium  isotopes  via 
a  3-step  ionization  process,  using  a  variety  of  input  parameters.2  As  seen  in  Figure  12,  the 
model  has  been  revised  to  perform  these  calculations  for  isotopes  of  plutonium  and  other 
isotopes  of  interest  as  prioritized  in  [15]. 


3rd  Excited  State 


Casimir 
(odd  isotopes) 


Transition  C 


2nd  Excited  State 


J3  ' 


J2 


Casimir 
(odd  isotopes) 


Transition  B 


1st  Excited  State 

J2 '  r 


J1 


Casimir 
(odd  isotopes) 


Transition  A 


Initial  State 

ji  :fo 


Casimir 
(odd  isotopes) 


AJo 

Bfo" 


MHz 

MHz 


nm 


Centroid 
[750.234 

FWHM 

|lE-10  nm 


Al0 

B|T 


MHz 

MHz 


nm 


Centroid 
i  847.275 

FWHM 

|l.7E-6  nm 


A|o 

Bfo" 


MHz 

MHz 


Centroid 

j  420.764  nm 

FWHM 

1 1.6E-6  nm 


A]0 

Bo 


MHz 

MHz 


A[o 
B  jo 


MHz 

MHz 


nm 


Centroid 
[750.239 

FWHM 

I  IE-10  nm 


A[o 

BfiT 


MHz 

MHz 


Centroid 

j  847.262  nm 

FWHM 

fl.7E-6  nm 


A[o 
B  jo 


MHz 

MHz 


Centroid 

1 420.77  nm 

FWHM 

1 1.6E-6  nm 


A[o 

b  nr 


MHz 

MHz 


239  242 

Figure  12.  Atomic  parameter  inputs  into  the  model  for  Pu  and  Pu. 


2  Assumes  the  atoms  are  in  the  ground  state. 
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As  shown  in  Figure  13,  the  model  begins  with  the  calculation  of  each  cross 
section  and  laser  spectral  irradiance  for  each  transition  we  are  attempting  to  decipher. 
Once  calculated,  they  are  combined  considering  the  time-dependent  laser  amplitude.  This 
yields,  “A  time-dependent  transition  rate.  These  rates  are  used  to  define  parameters  for 
the  rate  equations,  which  describes  the  rate  change  in  population  of  each  atomic  level” 
[8].  Upon  determination  of  the  parameters,  they  are  integrated  over  a  period  of  time,  and 
a  cumulative  fraction  of  ions  is  produced.  This  process  can  be  repeated  for  the  given 
isotope  and  iterated  over  other  values. 


For  each  transition: 


Figure  13.  Rate  equation  model  diagram,  from  [8]. 
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A. 


CROSS  SECTIONS 


1.  Assumptions 

The  cross  section  for  absorption  can  be  given  in  [16]  by 

An  27 ts(co) 

°m1m2=Tr21gM=1^\D12\2 

Isselhardt  defines  the  variables  as,  “A0  is  the  mean  transition  wavelength,  F21  is  the  partial 
width  of  the  excited  state,  g(w )  is  the  normalized  lineshape  of  the  transition,  and  \D12\2 
is  the  squared  magnitude  of  the  dipole  matrix  element  for  each  transition”  [8].  This  dipole 
element  can  be  expressed  as  the  product  of  a  reduced  matrix  element  and  the  Clebsch- 
Gordan  coefficient,3  per  the  Wigner-Eckart  Theorem  [8] 

{Jz^z  I D  even  1/i.mr)  =  C{]1,]2>A],m1,m2,Am'){]2\\D  even  ll/l) 

The  variables  are  defined  as:  “ C(J1,J2,AJ,m1,m2,Am )  is  the  Clebsch-Gordan  coefficient 
and  (J2  WDevenWJ-t)  is  the  reduced  matrix  element”  [8].  The  electron  spin  angular 
momentum  is  not  affected  by  the  dipole  operator. 

The  even  isotopes  have  equal  squares  of  their  Clebsch-Gordan  coefficients  for 
each  transition  between  each  substrate.  In  odd  isotopes,  the  transitional  energies  will  not 
be  the  same.  Using  the  selection  rules  stated  in  Table  2,  we  can  use  Wigner-Eckart 
Theorem  and  state  it  as  described  in  [8] 

(F2,m2\Dodd\F1,m1)  =  C(F1,F2lAF,m1,m2,Am)(J2\\Dodd\\J1) 

</2l|£WllA>=  (J2\\d  evenWJl ) 

2.  Transitional  Cross  Sections  for  Even  and  Odd  Isotopes 

The  model  calculates  the  cross  section  as  a  product  of  amplitude  and  nonnalized 
lineshape.  This  is  given  in  [8]  as 

flu  (A)  =  oh  •  giy ) 


3  Per  [17],  “Clebsch-Gordan  coefficients  arise  when  two  angular  momenta  are  combined  into  a  total 
angular  momentum.  These  coefficients  occur  when  angular  momentum  from  one  system  is  found  to  have  a 
combination  of  angular  momenta  from  two  subsystems.” 
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A2  . 

where  “er0  =  T21  i  s  the  peak  cross  section  and  g  (y)  i  s  the  lineshape  for  the  transition” 

[8].  This  lineshape  describes  combining  the  Doppler  broadening  and  width  due  to 
spontaneous  decay  and  is  the  union  of  Lorentzian  with  a  Gaussian  Functions.  Therefore, 
the  cross  sectional  for  even  isotopes  can  be  given  as  [16] 

^en=Tr2l---[®(A)®X(A)] 

4  9  i 

In  odd  isotopes,  the  differences  in  resonant  wavelength  of  each  transition  must  be 
accounted  for.  Isselhardt  explains,  “The  cross  section  is  proportional  to  the  square  of  the 
transition  dipole  and  can  be  multiplied  by  peak  cross  section  by  the  squared  Clebsch- 
Gordan  coefficient,  as 

omim2W  =  —  r2i  •  C(F1,m1,F2,m2y  ■  g(X) 

where  C^F^m-i,  F2,m2)2  is  defined  as  the  squared  Clebsch-Gordan  coefficient.  Since  in 
an  odd  isotope,  the  cross  section  is  proportional  to  the  sum  of  the  squared  Clebsch- 
Gordan  coefficient  for  each  transition”  [8].  The  transitional  energy  differences  due  to 
differences  in  total  angular  momentum  must  be  taken  into  account.  Therefore  from  [8], 


° odd  K 
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B.  LASER  IRRADIANCE 

The  spectral  irradiance  of  the  laser  / ,{  X  )  is  a  product  of  photon  flux  (</>,)  and  the 
laser  lineshape  (/(  X  )) 

hW  =  01  •  m 

Hurst  [14]  gives  the  photon  flux  as 

< p^photons  ■  cm~ 2]  =  (5  •  109)  •  (^[nm])  •  (P;[/r/])  /  A[cm2] 
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and  defines  Pi  as  the  pulse  intensity,  A  as  the  area,  and  as  the  wavelength  of  the  laser 
beam.  These  parameters  can  be  varied. 

According  to  Isselhardt,  the  intensity  of  a  laser  beam  can  be  characterized  by 
“Using  average  total  beam  energy  and  spatial  distribution  in  terms  of  the  1  /  e2  width” 
[8],  For  normal  Gaussian  distribution,  this  width  is  equivalent  to  4a.  Therefore,  it  can  be 
said  half  of  the  population  is  within  ±0.674a  of  peak  distribution  [8],  This  area  is  defined 
as  “^500/0  =  n  •  (.  674)rw  •  (.  674)rv,  where  rwand  rv  are  horizontal  and  vertical  radii. 
Also,  average  elliptical  area  is  given  by  A4a  =  n  •  16  •  rH  •  rv”  [8],  Using  this  knowledge, 
model  predication  and  experimental  results  can  be  compared  using  the  modified  laser 
intensity  equation 

Pavg=Pr  0.5--^-  =17.58-^ 

This  takes  experimental  values  of  laser  intensity  and  cross  sectional  area  and  converts  it 
into  average  pulse  energy  ( PaVg )• 


During  usage,  lasers  have  been  known  to  exhibit  variations  in  mean  wavelength, 
power,  position,  and  timing  from  pulse  to  pulse  [16],  In  order  to  reproduce  the 
fluctuations  in  wavelength  cause,  the  model  calculates  a  distribution  of  ionization 
probabilities  over  a  designated  wavelength  range.  It  then  randomly  selects  samples  in  this 
distribution  100,000  times  using  a  nonnal  distribution. 


Isselhardt  explains,  “The  model  was  built  around  the  rate  of  change  in  population 
of  atomic  level”  [8].  It  calculates  the  population  density  of  four  states  as  they  are 
irradiated:  ground  level,  two  bound  excited  states,  and  finally  the  ionization  level  [16]. 
These  four  equations  are  given  in  [8]  as 
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dNion 

dt 


—  ^2ion^2  +  W2cN2 


The  variables  are  defined  in  [8]  as:  Nt  is  the  number  of  atoms  in  state  |  i ),  WLJ  is  the  rate 
of  transition  from  state  |  i )  to  state  |  j  ),  gL  is  a  degeneracy  factor  for  state  |  i ),  and  Tt  is 
the  average  lifetime  of  state  |  i  ).  ±  W2CN2  terms  are  added  to  account  for  the  additional 
ionization  pathways.  The  negative  tenn  accounts  for  equation  of  state  |  2  )  and  the 
position  tenn  accounts  for  the  |  ion)  state  [8].  Once  the  rates  WLj  are  defined  by  [8]  as 

Wij  =  f  f  °ijW  ■  Ji(A,  t)  dA  dt  =  ^  ^  Oij{ 4)  •  7j(4,  t)AXAt 

l  k 

the  rate  equations  can  be  solved.  cri;(A)  is  the  cross  section  between  |  i  )  and  j  ),  and 
7;(/l,  t)  is  spectral  inadiance  of  the  laser  used  to  cause  the  excitation.  If  a  continuous  laser 
was  used,  we  would  integrate.  Since  a  pulsed  laser  is  used,  we  use  the  summation. 

This  rate  depends  on  a  time  variance  in  the  amplitude  of  the  laser  pulse.  However, 
this  time  dependence  is  independent  of  wavelength  and  is  expressed  in  [8]  as 

1 Wij{t)  =  Wij  ■  TL(t) 

where  the  variables  expression  is  defined  as  'Tj(t)  being  the  time  distribution  of  the  laser 
pulse  [8].  The  lineshape  is  assumed  to  be  Gaussian  in  nature  and  is  described  as  [16] 


Tt(t)  = 


fi-r0)2 

2a2 


where  Isselhardt  defines  the  variables  as:  a  is  the  standard  deviation  and  T{)  is  the  time 
taken  at  the  peak  of  the  pulse. 
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IV.  EXPERIMENTAL  AND  MODEL  RESULTS 


In  the  effort  to  use  the  model  to  make  comparisons  with  experimental  data,  some 
troubleshooting  had  to  occur  in  order  to  obtain  the  correct  reference  data.  A  trial  and  error 
process  was  performed  on  both  the  first  and  second  resonance  laser,  but  only  the  second 
is  shown  in  the  description  below.  Section  A  documents  this  troubleshooting  process. 
Table  3  provides  a  list  of  all  parameters  used. 

A.  TROUBLESHOOTING 

The  troubleshooting  effort  consisted  of  six  tests.  Three  of  them  are  documented 
here.  The  tests  consisted  of  a  trial  by  error  method  of  determining  the  most  suitable  value 
of  a  parameter.  They  were  based  on  previously  known  data  for  both  the  first  and  second 
resonant  laser.  This  process  was  aimed  at  determining  values  for  integration  limits  and 
laser  power,  ensuring  the  end  plot  was  similar  to  Figure  16. 

Figure  14  presents  test  1,  which  is  a  modeled  version  of  the  detennination  of  the 

242  239 

Pu/  Pu  ratio  vs.  wavelength  of  the  second  resonance  laser.  As  depicted  below,  it  can 
be  seen  that  as  the  laser  bandwidth  was  broadened,  the  relative  lineshape  was  also 
broadened.  As  can  be  seen  in  the  figure,  the  results  for  a  bandwidth  of  0.00 lnm  are  not 
visible  on  the  plot,  which  led  to  further  investigation  to  detennine  the  reason  for  this. 
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Figure  14.  Test  1.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  second 

resonance  laser  at  full  power. 


Figure  15  shows  the  results  for  test  2,  another  modeled  version  of  Pu/  Pu  ratio 
vs.  wavelength  of  the  second  resonance  laser.  Unlike  test  1,  test  2  assumed  a  reduced 
laser  power.  The  changes  from  the  test  1  results  clearly  illustrate  how  laser  power  effects 
the  calculation  of  isotope  ratio.  The  lineshape  for  bandwidth  O.OOlnm  is  now  visible  in 
this  plot.  However,  the  peaks  for  all  bandwidths  centered  around  847.195nm-847.22nm 
and  847.32-847.35,  appear  anomalous  and  should  not  be  included. 
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Figure  15.  Test  2.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  second 
resonance  laser  with  laser  power  reduced  by  a  factor  of  ten. 


Figure  16  presents  the  results  of  test  3.  Laser  power  was  maintained  at  its  reduced 
level,  but  the  integration  limits  inside  the  model  were  expanded.  As  illustrated,  the 
anomalous  peaks  at  the  end  and  beginning  of  the  each  bandwidth  have  disappeared. 
These  line  shapes  are  more  representative  of  what  is  expected.  The  values  of  Table  3,  test 
3  were  used  in  the  remaining  tests,  except  when  noted. 

This  entire  process  demonstrates  the  ability  of  the  model  to  easily  troubleshoot 
and  update  fields  for  experimental  data,  without  having  to  unnecessarily  utilize  precious 
resources.  Figure  12  displays  an  example  of  the  fields  in  which  the  ability  to  input  a 
desired  value  is  available. 
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Figure  16.  Test  3.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  second 
resonance  laser  with  laser  power  reduced  by  a  factor  of  ten  and 
expanded  integration  limits. 
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Table  3.  Parameters  used  in  troubleshooting.  Laser  power  is  reduced  by  a 
factor  of  10  in  test  2,  while  maintaining  the  same  integration  limits. 
Test  3  maintains  the  reduced  laser  power,  but  expands  the  integration 
limits.  Data  was  drawn  from  [18]. 


Test  1 

Parameter 

Laser  Power  [mW] 

Integration 
Limits  [nm] 

First 

Second 

Third 

Resonant 

Resonant 

Resonant 

1000 

656 

200 

847.15- 

847.38 

Test  2 


Parameter 

Laser  Power  [mW] 

Integration 
Limits  [nm] 

First 

Second 

Third 

Resonant 

Resonant 

Resonant 

100 

65.6 

20 

847.15- 

847.38 

Test  3 


Parameter 

Laser  Power  [mW] 

Integration 
Limits  [nm] 

First 

Resonant 

Second 

Resonant 

Third 

Resonant 

100 

65.6 

20 

847-847.75 

B.  EXPERIMENTAL  AND  MODEL  COMPARISONS 

Figure  17  shows  previous  experimental  results  from  [18]  for  three  values  of  laser 
bandwidth  (0.00 lnm,  0.005nm,  and  0.015nm).  Figure  18  shows  modeled  results  for  the 
same  three  bandwidths.  In  larger  bandwidths,  the  isotopes  are  saturated  in  between  the 

peaks.  This  is  not  the  case  for  the  smaller  bandwidths.  It  is  also  clear  to  see  that  for  the 
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0.015nm  bandwidth,  the  peaks  are  more  symmetric  than  the  0.00 lnm  bandwidth. 
According  to  Isselhardt,  this  is  caused  by  the  fact  that  “the  tails  of  the  isotope  not 
saturated  in  that  wavelength  region  of  the  peak  become  less  sensitive  to  wavelength”  [8]. 
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Figure  17.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  first  resonance  laser 
using  experimental  laser  parameters,  from  [18]. 
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Figure  18.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  first  resonance  laser 

using  modeling  parameters,  from  [19]. 
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Figure  19  shows  additional  experimental  results  from  [18]  for  three  values  of  laser 
bandwidth  (0.00 lnm,  0.005nm,  and  0.015nm).  Figure  20  shows  modeled  results  for  the 
same  three  bandwidths.  This  indicates  an  isotope  ratio  where  each  resonance  is  on  either 
side  of  the  peak  and  they  are  about  equal  in  amplitude  where  they  cross  the  y-axis  value 
of  one  at  approximately  847.529nm  on  Figure  19  and  847.27nm  on  Figure  20.  However, 
Figure  20  differs  from  19  in  that  each  peak’s  amplitude  is  approximately  a  factor  of  ten 
larger.  It  also  seems  that  the  peaks  for  each  bandwidth  are  reversed. 
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Figure  19.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  second  resonance 
laser  using  experimental  laser  parameters,  from  [18]. 
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Figure  20.  Isotope  ratio  of  Pu/  Pu  vs.  wavelength  of  the  second  resonance 

laser  using  modeled  parameters. 


C.  SUMMARY 

The  research  in  this  thesis  included  a  series  of  troubleshooting  tests  to  provide 
calibration  and  understanding  of  the  impact  of  parameter  choices  to  improve  the  fidelity 
of  the  modeling  application.  Then  it  involved  exercise  of  the  model  to  compare  outputs  to 
known  experimental  data  to  gain  important  insights  into  the  application  of  the  model  to 
problems  of  significance  to  the  RIMS  nuclear  forensics  assessment. 
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V.  CONCLUSION 


It  is  believed  that  9  countries  have  nuclear  weapons,  with  others  being  considered 
primed  for  development  of  such  a  capability  [20].  Each  of  these  nuclear  weapons 
programs  has  distinguishable  characteristics  defined  by  their  individual  processes.  If  a 
nuclear  attack  were  to  occur,  a  most  suitable  and  expeditious  method  for  detennining  the 
origin  of  the  device  or  its  nuclear  materials  is  the  RIMS  technology.  There  are  limited 
resources  for  experimental  determination  of  the  data  needed  for  RIMS  analysis.  With  that 
in  mind,  the  rate  model  equation  utilized  in  this  research  is  established  as  a  necessity.  We 
can  use  this  model  to  assist  in  predicting  ionization  probabilities  for  many  elements  and 
their  isotopes,  without  depleting  the  reserve  of  resources. 

A.  SUMMARY 

The  research  started  with  a  troubleshooting  process.  It  was  known  how  the  plot’s 
lineshape  should  look.  Various  values  of  laser  centroid,  power,  and  integration  limits 
were  chosen  until  the  known  plot  from  experimental  data  was  replicated  by  modeling 
results.  Such  a  plot  is  displayed  in  Figure  16,  with  the  values  of  each  parameter  that 
resulted  from  the  troubleshooting  process  in  Table  3.  Once  these  values  were  detennined, 
the  comparing  and  contrasting  of  experimental  and  model  parameters  commenced. 

As  seen  in  Figures  17  and  18,  the  rate  equation  model  was  able  to  sufficiently 
replicate  experimental  data.  It  consistently  displayed  similar  characteristics  to  what  has 
been  seen  experimentally  for  the  isotope  ratios  (  Pu/  Pu)  when  excited  by  the  first 
laser  wavelength.  Displayed  in  Figures  19  and  20  are  the  plots  of  the  modeled  and 
experimental  laser  parameters  for  the  isotope  ratio  versus  wavelength  of  the  second  laser. 
It  can  be  clearly  seen  that  a  disagreement  between  model  and  experimental  data  is 
present.  This  indicates  that  more  research  is  necessary  in  order  to  reproduce  acceptable 
isotope  ratios  for  the  second  resonance  laser  wavelengths.  However,  this  is  not 
considered  to  detract  from  the  present  research;  it  implies  that  the  best  way  to  interpret 
and  utilize  the  model  predictions  is  to  compare  the  results  to  more  precise  experimental 
data  and  continue  the  process  of  fine  tuning  the  model. 
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B. 


FUTURE  WORK 


In  Figure  21,  the  period  table’s  actinide  series  elements  have  been  highlighted. 
Per  [15],  americium  (Am)  and  neptunium  (Np)  are  the  next  elements  to  be  researched. 
Both  elements  “are  primarily  formed  from  the  decay  of  Pu  and  U”  [15].  Understanding 
these  elements  and  how  their  isotopes  are  formed  will  assist  in  detennining  the  fuel 
isotopes  and  where  the  fuel  was  manufactured.  Np  ([21],  [22])  and  Am  ( [23]— [25])  have 
already  had  preliminary  studies  to  determine  some  of  the  parameters.  However,  in  order 
to  fully  experiment,  further  research  is  required  to  gain  key  parameters. 

Work  is  also  being  conducted  to  substitute  one  of  the  pulsed  lasers  with  a 
continuous  wave  (CW)  laser.  It  would  be  beneficial  to  re-run  some  of  the  experiments  on 
U  and  Pu  with  the  new  laser  implemented.  Running  the  model  based  on  the  parameters  of 
the  CW  laser  would  be  valuable  in  validating  and  understanding  the  experiments. 


PERIODIC  TABLE  OF  THE  ELEMENTS 


H 

He 

1.0079 

4.0026 

3 

4 

5 

6 

7 

8 
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Figure  21.  Periodic  table  of  elements,  from  [26].  Highlighted  in  red  are  the 
isotopes  in  the  actinide  series  in  which  we  are  most  interested. 
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